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1 Introduction 

Schwinger's original nonet mass formula [IJ (here the symbol for the meson stands 
either for its mass or mass squared), 

(AK - 3r] - tt) (377' + 7T - AK) = 8(K - tt) 2 , (1) 

and the Sakurai mass- mixing angle formula 0, 

2 AK - 3r) - 71 
tan 9 = 7F' 2 

both relate the masses of the isovector (tt), isodoublet (K) and isoscalar mostly octet 
(77) and mostly singlet (77 ) states of a meson nonet, and the nonet mixing angle (9). 
We alert the reader that, although we use notation suggestive of masses below, each 
formula is to be reprised in terms of mass or mass squared values in this introductory 
section. 

The relations (1) and (2) are usually derived in the following way: For a meson 
nonet, the isoscalar octet-singlet mass matrix, 

is diagonalized by the masses of the physical 77 and 77 states: 

/ cos 9 sin 9 \ I rj \ / cos 9 — sin 9 \ 



I — sin 9 cos 9 J \ rj J I sin cos 1 

/ cos 2 6* 77 + sin 2 6* 77' sin 9 cos ^ (77' — 77) \ 
I sin 9 cos # (77 — rj) sin 2 # 77 + cos 2 9 r] j 



(4) 



where 9 is the nonet mixing angle, which is determined by comparing the correspond- 
ing quadrants of the matrices (3) and (4), by any of the three following relations: 

, 2 a M 88 -r] 

tan 9 = 7^4? (5) 

tan 2 £ = ^-^, (6) 
M 99 - 77 V ; 

2M«q 

sin2^ = 4^- (7) 
77—77 

It is easily seen that Eqs. (5) and (6) are identical, since, due to the trace invariance 
oiM, 

77 + 77' = M 88 + M 99 , (8) 

and therefore, M 88 — rj = rj' — M 99 , and 77' — M 88 = M 99 — 77. Eliminating 9 from 
(5), (7), or (6), (7), with the help of sin 29 = 2 tan 9/ (1 +tan 2 9), leads, respectively, to 

(77-M 88 )(M 88 -t/) = M 2 9 , (9) 
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(M 99 - V )( V -M 99 ) = M 2 9 , (10) 

which again are identical, through (8). 

We note that the "ideal" structure of a meson nonet, 

(rj = 2K — it, i] = 7r), or (rj = it, rj = 2K — tt), (11) 

and the corresponding pure qq physical state valence flavor wavefunctions, 



uu + dd 

-ss, 0J n > = 

given by the mixing 



LUrj = —ss, uj i = -i= — = nn, or u v = nn, to > = ss, (12) 

V2 



u n = u 8 cos9 — u 9 sm9, 

u> = u)g sin 9 -\- ujg cos 9 (13) 
with the "ideal" nonet mixing angle, 

9 = arctan = 35.3°, or 9 = - arctan V2 S -54.7°, (14) 

and the definitions 

uu + dd — 2ss uu + dd + ss . . 

-s = , C = , (15) 

are the unique solution to Eqs. (9), (10) under the quark model inspired conditions 

M 88 = l^, Af 99 = ^, M m = - 2 -^(K-*), (16) 

where the first of the three relations in (16) is the standard Gell-Mann-Okubo mass 
formula ||. 

For all well established meson nonets, except the pseudoscalar (and, we expect, 
scalar) one(s), both linear and quadratic versions of Eqs. (5)-(10) are in good agree- 
ment with experiment. For example, for vector mesons, if one assumes the validity of 
the Gell-Mann-Okubo formula u)$ = (4K* — p)/3, then one obtains from Eq. (9) with 
the measured meson masses H, M m = —0.209 ± 0.001 GeV 2 in the quadratic case, 
and —0.113 ± 0.001 GeV in the linear case. Note that this is entirely consistent with 
-0.196 ± 0.005 GeV 2 and -0.118 ± 0.003 GeV, respectively, which follow from the 
third element of (16). For tensor mesons, a similar comparison gives —0.305 ± 0.020 
GeV 2 vs. -0.287 ± 0.015 GeV 2 , and -0.105 ± 0.008 GeV vs. -0.104 ± 0.005 GeV, 
respectively. 

However, for the pseudoscalar nonet, one obtains from Eq. (5) with meson masses 
squared, 9 w —11°, in sharp disagreement with experiment, which favors the 77-77 
mixing angle in the vicinity of —20° [|], [5], Although using linear meson masses 
in Eq. (5) does give 9 w —24°, in better agreement with data than its mass-squared 
counterpart, the value of M 89 , as given by (9), is now —0.165 ± 0.004 GeV, vs. 
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—2y/2/3 (K — 7r) = —0.338 ± 0.004 GeV. This emphasizes that neither Schwinger's 
nonet mass formula nor the mass- mixing angle relations (including Sakurai's) (5)- (7) 
hold for the pseudoscalar nonet. 

It is well known, however, that the pseudoscalar (and, probably, scalar) mass 
spectrum does not follow the "ideal" structure, Eq. (16), since the mass of the 
pseudoscalar isoscalar singlet state is shifted up from its "ideal" value of (2K + n)/3, 
presumably by the instanton-induced 't Hooft interaction J7| which breaks axial U(l) 
symmetry However, the use of M 99 = (2K+ir)/3+A, A ^ 0, in Eqs. (5)-(8) 

will again lead to Schwinger's formula (9), which does not hold for the pseudoscalar 
mesons, as just demonstrated. [In fact, the structure of this formula does not depend 
at all on M 99 , as seen in (9).] Therefore, instanton, as well as any other effects which 
may shift the mass of the pseudoscalar isoscalar singlet state, cannot constitute the 
explanation of the failure of Schwinger's quartic mass and the Sakurai mass-mixing 
angle formulae for the pseudoscalar nonet. We believe, however, that the following 
analysis can resolve this problem. 



2 Pseudoscalar meson mass squared matrix 

It is known that the observed mass splitting among the pseudoscalar nonet may be 
induced (in terms of the 1/N C expansion) by the following symmetry breaking terms 



= L (bTiM(u + U^) +4t [TrflnCZ-lnC/ 1 
with M being the quark mass matrix, 
M = diag (m u , m d , m S/ 



(17) 



m s diag (x,y, 1), x = — ,y= — , 



N c the number of colors, / the pseudoscalar decay constant in the limit of exact nonet 
symmetry, and B,e = const. This symmetry breaking Lagrangian term is to be added 
to the U(3)lxU(3)k invariant non-linear Lagrangian 



f 2 

L (o) = J Tr (q^uquu^ 



(19) 



with 



U 



0,1,. ..,8, 



exp(z7r//), 7r = A a 7r a , a 
which incorporates the constraints of current algebra for the light pseudoscalars ir a 

As pointed out in ref. [12|, chiral corrections can be important, the kaon mass 
being half the typical 1 GeV chiral symmetry breaking scale. Such large corrections 
are clearly required from the study of the octet-singlet mass squared matrix M 2 . In 
the isospin limit x = y one has [13|] (with m n = {m u + m d )/2) 



M 



B 



|(2m s 



m. 



^(m n - m s ) 



^(m n -m s ) |(m s + 2m n ) + ^ 



(20) 
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which, on the most naive level, through the schematic Gell-Mann-Oakes-Renner re- 



lations (to first order in chiral symmetry breaking) [14] 



7T 

K 2 

v 2 



2B m n , 

B (m s + m n ), 
2 

-B {2m s + m n ), 



(21) 



which we discuss in more detail below, reduces to 



M 



AK 2 - it 2 

-2V2 (K 2 - 7T 2 ) 



-2V2 (K 2 - 7T 2 ) 

2K 2 + n 2 + 3A 



A 



BN r 



(22) 



Also, the Nambu-Jona-Lasinio model with the instanton-induced 't Hooft inter- 
action, initiated by Hatsuda and Kunihiro |T3| and Bernard et al. 



16 , and then 



extensively studied by Dmitrasinovic || [TCfl, provides shifts of both the pseudoscalar 
and scalar isoscalar singlet masses by the same amount, but in opposite directions 
(viz., the pseudoscalar isoscalar singlet mass is increased, while the scalar one is de- 
creased) . Thus, for the pseudoscalar mesons, rather than using the model-dependent 
A, as defined in (22), we introduce the quantity A below, which may be considered 
as the sum of all possible contributions to the shift of the isoscalar singlet mass (from 
instanton effects, l/iV c -expansion diagrams, gluon annihilation diagrams, etc.). 

Here, however, we suggest that the actual form of the mass squared matrix for 
the pseudoscalar mesons (and the corresponding symmetry breaking terms in (17)) 
is as follows: 



f 2 M 2 



±PkK 2 - f 2 n 2 
-2V2 (f 2 K K 2 - /> 2 ) 



-2V2 (f 2 K K 2 - f 2 n 2 ) 
VlK 2 + f 2 ^ 2 + 3f 2 A 



(23) 



where /'s are the pseudoscalar decay constants defined below. 

Indeed, the form of such a mass squared matrix is determined by the form of 
Gell-Mann-Okubo type relations among the masses of the isovector, isodoublet, and 
isoscalar octet and singlet states (which in our case are Eqs. (26)- (28) below), since 
this matrix must be equivalent to that of the form (3), which in the case we are 
considering is 

( fivis hhrfw \ 



f 2 M 2 



(24) 



and is diagonalized by the physical 77 and 77 meson masses and decay constants: 

\ 



f 2 M 2 



/„V 



(25) 



The equivalence of the matrices (23) and (24) is guaranteed by the validity of the 
following relations: 



2„2 



1 r 



m u (uu) + md{dd) + Am s (ss) 



m* 2 - fw 



(26) 
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a/2 r - i 2\/2 

fsfdVsa = — 5- m u {uu) + m d (dd) - 2m s (ss) 



fdVgg = H A - - m u (uu) + m d (dd) + m s (ss) = f g A + 



O 2 - f K K 2 

viK 2 + ry 



with 



K 2 



(X±) 2 + (K°) 



(27) 
(28) 

(29) 



as suggested by Dmitrasinovic [ TO | , on the basis of the (precise) Gell-Mann-Oakes- 
Renner formulae which relate the pseudoscalar masses and decay constants to the 
quark masses and condensates ||l4j| : 



T *r 7T 

J 7T 

2 



— m u (uu) + m ( i(dd) 

— [m u (uu) + m s (ss)] 

— md(dd) + m s (ss) 



(30) 
(31) 
(32) 



(We ignore Dashen's theorem violating effects JT7|, and only approximately take into 
account isospin violating effects via (29), as we are not concerned here with accuracies 
better than 1%.) 

Note that in the limit of exact nonet symmetry, 



U = Ik = fss = fm = f, (uu) = (dd) = (ss) = (qq), 



(33) 



one has the mass squared matrix (20) with B = —(qq)/f 2 , which further reduces to 
(22). The real world, and the mass squared matrix (23) associated with it, however 
(as we shall see), corresponds to the situation when the first set of the relations (33) 
is broken, but the second one remains (approximately) satisfied, i.e., the amount of 
SU(3) flavor symmetry breaking in terms of the quark condensates is much smaller 
than that in terms of the pseudoscalar decay constants. 



3 Modified Gell-Mann— Okubo mass formula 

Here, we shall only explicitly demonstrate the validity of the modified Gell-Mann- 
Okubo formula (26) (the remaining relations (27) and (28) may be checked in a similar 
way). 

First, this relation, as well as the mass squared matrix (23), (24), may be obtained 
through relating the vacuum expectation values of the equal time axial divergence- 
axial current commutators (which are "sigma" commutators in Eq. (34), (35) below) 
to integrals over the pseudoscalar meson spectrum, as done by Gensini |TH] following 



Gatto et al. [19]. For the symmetry realized through a set of massless Goldstone 
pseudoscalar mesons, only the pole terms survive in the first order of the symmetry 
breaking and are expected to dominate over the continuum, which contributes only 



at the second-order level. We therefore have the identities [ITU 



<M0)> = faf b M 2 ab + / d/i Pab (fi), (34) 
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which define the current quark mass contributions to the pseudoscalar masses. These 
may be decomposed as 



M 



2 

ab n n 

Jajb 



where M^ m is the long-range electromagnetic contribution, Mg is the gluon term 
which stands for both perturbative two-gluon annihilation and nonperturbative in- 
stanton effects. (Both are responsible for the shift of the isoscalar singlet mass prior 
to its mixing with the isoscalar octet, leading to the physical rj and rj masses.) The 
latter must therefore act only on the singlet-singlet state matrix element: 

(M*) ab = A8 a9 8 b9 . (36) 

Finally, M| o t stands for the contribution of the higher order terms. 

Mass terms, which are the only explicit symmetry-breaking terms, have the general 
formQ 

L m (x) = £ S {x) + £ 3 S 3 (x) + . . . + £ A r2_ 1 S , A f2_ 1 (x), (37) 

where S { = ip(x) \/2 ip(x), A = (2/iV) 1/2 I, and X M ^ _i are the Gell-Mann-type 
matrices which constitute the SU(M) basis, M < N. Restricting the indices to a 
subgroup SU(iV ), N < N, does not introduce higher mass quark fields anywhere but 
in the singlet-singlet part, where they can be absorbed into the gluon term within 
the definition of higher-order terms given above. This restriction is meaningful only 
if the mixing of low mass ipip pahs with higher mass ones is small enough. This is 
precisely what preserves the approximate SU(3) symmetry of hadronic interactions, 
independent of the total number of flavors. 

In terms of the current quark masses we also have 



, N .-1/2N 

£ o = [yj z2 m ^ ( 38 ) 

/M-l \ 

e M 2_ 1 = [M(M-l)]- 1/2 f J2m { - (M-l)m M j. (39) 

Assuming the SU(3)-invariant vacuum, fixing M — 3, and neglecting long-range elec- 
tromagnetic effects, one obtains in the lowest order of the symmetry breaking 



(**) ! 



P 

— B (m u + m d ), 



JK 



P 

Pk 



^ B (m d + m s ) (40) 



1 We use the notations of ref. mm 
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in the charged sector, and 



^33 



^38 



^39 



^88 



^89 



^99 



P 

— B (m u + m d ), 

P ft 111 **- m d 



V3 ' 



P_ 

P 4m s + m u + m d 

Tl B 3 • 

/8/9 3 

P 2 

~2 B - (m s + m u + m d ) + A 

J9 



(41) 



in the neutral Y = 1% = sector. Eliminating the quark masses from the above two 
sets of relations, one obtains 



r-2 2 

r2 2 
i8 %8 

p2 2 

h V99 



ft M 

\[2fl 
V2 



K ± ) + 2{K' 
Pk + (*°) 



PA* 



1 

3 



Pk^Y + (*°) 2 ) + P 



7T 



+ PA 



(42) 



which corresponds to Eqs. (26)- (28). Using the definition (20), one can also obtain the 
mass squared matrix (23), (24) from Eqs. (40)- (42). (We have, consistently, ignored 
small 7T° — i] and rj [7738,7739] mixing effects of electromagnetic and isospin breaking 
origin). 

Independently, the relation (26) may be obtained from the following expressions 
for the pseudoscalar meson masses calculated by Li [2D| in the chiral effective field 
theory of mesons (fi, A = const): 



7T 



(*•)= 



16iV cy u 3 /A 2 \ . . 

U 7r )2 f 2 I ln — - 7 + 1 I ( m u + m d ) 

(47r)2/ | \ ^—-l+l)(rn u + m s 



A 2 
,2 



vis 



16N c [m 
W?P 



K \ A* 

3 ' A 2 \ , 

ln — - 7 + 1 [m d + m s 
k \ A 4 / 

3 / A2 

In — - 7 + 1 



m u + m d + 4m s 



(43) 



N.S 
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Additionally, even at lowest order, 0(p 2 ), in generalized chiral perturbation theory, 
the inclusion of decay constants is required to form a relation among the octet pseu- 
doscalar mesons ||2T|. Finally, this formula may be also obtained in standard chiral 
perturbation theory f22|| , as follows. 

Standard chiral perturbation theory leads to the following expressions for the 



pseudoscalar meson masses and decay constants [22 



7r 2 = 2m n B f 1 + ii n - -jj 8 + 2m n K 3 + K 4 
( 2 

K 2 = (m n + m s )B f 1 + -/i 8 + (m n + m s )K 3 + K± 

2 / 4 2 \ 

^88 = g( m ™ + 2m s )B M + 2fj, K - -fi 8 + ~{m n + 2m s )K 3 + KA 

+ 2m n B + ~/j, K + ~/j,^J + K 5 , 

U = 7(1- 2/i,r - fix + 2m n K 6 + K 7 ) , 

_ / 3 3 3 \ 
fx = f ( 1 - ~ -I^k ~ -/is + {m n + m s )K 6 + K 7 j , 

f 8 = f(l-3tiK + l(m n + 2m s )K 6 + K r y (44) 

where /z's are chiral logarithms, and the constants K{ are proper combinations of 
the low energy coupling constants Lj. It then follows from these relations that the 



standard Gell-Mann-Okubo formula is broken in first nonleading order |22 



A GMO = 4iT 2 - 3r/ 2 8 - tx 2 = -2 (AK 2 !i K - 3 % 2 8/ u 8 - ix 2 ^) + • • ■ 



AB 



m n {n* + /i 8 - 2\i K ) +m s (2fi 8 - 2/i K ) +..., (45) 



where . . . stands for the higher order terms. 

However, the modified Gell-Mann-Okubo formula remains valid in this order, and 
is violated only by second order SU(3)-flavor breaking effects: 

a gmo = p (^IrK 2 - 3/ 8 2 r^ 8 - / 2 7T 2 ) = AB (m s - m n ) (p K + ~// 8 - -y^j + ... , 

(46) 

since the second factor on the r.h.s. of (46) must vanish in the SU(3)-flavor limit. 
Thus, the above analyses show, in an almost model independent way, that the modi- 
fied Gell-Mann-Okubo formula (26) is the only valid relation among the octet pseu- 
doscalar mesons,^ and therefore, the form of the mass squared matrix (23) is com- 
pletely justified. f\ 

2 A search for a relation of a more general form, Af^K 2 — 3/g?7| 8 + /"7r 2 , which would hold in 
the first nonleading order of standard chiral perturbation theory, results in a = 2. 

3 It has been suggested in the literature that the pseudoscalar decay constants should enter rela- 
tions like (26)-(28) in the first rather than second power J23|. As discussed above, such relations are 
expected to be less accurate than ours, according to chiral perturbation theory. 
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4 The Schwinger and Sakurai formulae reexam- 
ined 



Starting with the mass squared matrix (23), the considerations which lead to Eqs. 
(5)-(10) above, will now lead, through (26)-(28), to the following two relations, 



(4/^ 2 - 3/,y - jy) (s/ 2 , ^ + fy - 4 ^ 2 ) = 8 (/|^ 2 - ty) 2 , (48) 

which we refer to as "the Schwinger nonet mass, and Sakurai mass-mixing angle 
(respectively), formulae reexamined". 

In contrast to f n and fx, the values of which are well established experimentally 



y/2f K = 159.8 ± 1.6 MeV, y/2f n = 130.7 ± 0.3 MeV, ±f- = 1.22 ± 0.01 (49) 

the values of f v , ft and 9 are known rather poorly. We now wish to calculate the 
values of /,,, / / and 9, using the relations (47), (48), and compare the results with 
available experimental data. It is obvious that the two relations are not enough for 
determining the three unknowns. Note that the additional relation, independent of 
(47), (48) (the trace condition for (23), (25)), 

fW + fy = VI K 2 + fiA, (50) 

introduces an additional unknown, A. We therefore develop another independent 
relation among f v , f t and 9, as follows. 

The light pseudoscalar decay constants are defined by the matrix elements 

(0|$(0)7VyV>(0)|P(p)> = iS jP fpp", (51) 

where ip = (u, d, s) is the fundamental representation of SU(3)/, and P = r/ 88 , 7799) 
means the corresponding SU(3) / indices 3,8,9, thus picking out the diagonal (j = 3, 8) 

SU(3)/ Gell-Mann matrices A J , and A 9 = v/2/3 I. The neutral pseudoscalar wave 
functions, P, may be expressed in terms of the quark basis states qq: 

l^=E^k9)=E<l9^ ? = M,«, (52) 

g q 

where for P = vr°, c 3 u = l/y/2 = XlJy/2 = -X 3 22 /V2 = -c% c 3 s = 0, for P = r/ 88 , 
ct = c 8 d = 1/VE = X 8 U /V2 = XIJV2, c* = -2/V6 = AI3/V2, and for P = Vm 

cl = cl = ci = 1/V3 = (xyv2) qq . 
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The pseudoscalar decay constants denned in (51) can now be expressed as 

/p = E ^f-f*, (53) 

where we have introduced the auxiliary decay constants f qq - defined as the decay 
constants of the qq pseudoscalar bound states having the masses M(qq). In the isospin 
limit, f uU = fdd = fud = f-K = /tt+- Using this approximation, and evaluating the 
appropriate matrix elements leads to the following relations: 

,'cosfl - V2sme\ 2 , [sm8 + y/2cos6\ 2 , tr .. 
fv = ( at ) f* + I J fss, (54) 

2 2 

i ■ v - i . (cos9-V2^n9\ 
fv' = I TR I f« + [ ^ ) ( 55 ) 

Now we have four equations, (47), (48), (54), (55), which allow us to determine three 
unknowns, /,,, / / , 9, as well as the additional quantity introduced, namely, f sS . The 
solution to these four equations is 




fv 

u 

U 

fss 

9 = (-21.4 ± 1.0)°. (59) 



1.085 ±0.025, (56) 
1.195 ±0.035, (57) 
1.280 ±0.060, (58) 



[The 7r and K electromagnetic mass differences, and the uncertainties in the values 
of fa and fx, (see (49)), are taken as a measure of the uncertainties of the results.] 

Before comparing the solution obtained with experiment, let us also calculate the 
values of / 8 and f 9 which are obtained from (54), (55) in the no-mixing case (9 = 0): 

fs = \U + \U (60) 
h = \U + T^fss- (61) 



Therefore, as follows from (58), (59), 

u 

h 



1.185 ±0.040, (62) 
1.095 ±0.020. (63) 



The 77-77' mixing angle, as given in (59), is in agreement with most of experimental 
data which concentrate around —20° [§, || H- Also, the values for fs/ fn, hi fir and 
9 are consistent with those suggested in the literature, as we show in Table I. 
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Ref. 




fo/f* 


9, deg. 


This work 


1.185 ±0.040 


1.095 ±0.020 


-21.4 ± 1.0 


[6] 


1.11 ±0.06 


1.10 ±0.02 


-16.4 ± 1.2 


[25,26] 


1.25 


1.04 ±0.04 


-23 ±3 


[27] 


1.33 ±0.02 


1.05 ±0.04 


-22 ±3 


[28] 


1.12 ±0.14 


1.04 ±0.08 


-18.9 ±2.0 


[29] 


1.38 ±0.22 


1.06 ±0.03 


-22.0 ±3.3 


[30] 


1.254 


1.127 


-19.3 



Table I. Comparison of the values for fs/ fn, fvl ' fn and 9, calculated in the paper, 
with the results of the papers referenced. 



Note that (62), (63) are almost identical to (57), (56), respectively. This suggests 
that the quark content of the states must be the same. However, the relative phases, 
of course, cannot be. As we have suggested elsewhere this conundrum may be 
resolved by identifying the wavefunctions in (15) with 77 and 77 respectively, but with 
the signs of the ss terms of each reversed, viz., 

uu + dd — ss 1 uu + dd + 2ss 

" M — 73— • " M — — ' 



5 Comparison with data 

We now wish to compare the values obtained above for the ratios fs/ fn, f%l fn, and 
for the 77-77 mixing angle with available experimental data. We shall first consider 
in more detail the well-known 7r°, 77, 77 ^77 decays, for which experimental data 
are more complete than those for other processes involving light neutral pseudoscalar 
mesons, and then briefly mention the 77, 77 — » 7r + 7r~7, and J/ip — > 777, r\ 7 decays. 



5.1 P° — > 77 decays 

In the case of 7r° — > 77, the anomalous Wess-Zumino-Witten chiral lagrangian predicts 



31 



Aro^n = ~ ^e^e 2u k la k 2/3 , (64) 

37T/ 

where a = 1/137.036 is the fine structure constant, and e i: k\ are the polarization and 
momenta of the outgoing photons. A leading-log calculation of the chiral corrections 
reveals that the dominant effect is simply to replace / by the physical value f n 



The resulting amplitude is guaranteed by general theorems to remain unchanged in 



higher chiral orders p3| . One then finds that the predicted amplitude, as extracted 
from Eq. (71) below, with the experimentally measured width [|J r(7r° — > 77) = 
(7.7±0.6)eV, 

F^ 77 (0) = ^ = 0.025 GeV-\ (65) 
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is in excellent agreement with experiment 

F^ 77 (0) = (0.025 ± 0.001) GeV -1 , 



(66) 



thus providing the confidence that one may analyze the corresponding 77, 77' decays 
with a similar precision. 

In the case of the 77, 77 — > 77 decays, one should include both the 77-77 mixing 
and the renormalization of the octet-singlet couplings, which leads to the predicted 
amplitudes 



1/77(0) 



- f^cos^-2v / 2^ sin9 

n \J8 J 9 j 

- (^sin^ + 2 v / 2^cos^ 

fir \J8 J 9 j 



The values of these amplitudes, as extracted from data, are [|8|] 

F W7 (0) = 0.024 ± 0.001 GeV" 1 , 
F r? ' ry (0) = 0.031 ± 0.001 GeV" 1 . 

Calculation with the help of Eqs. (40), (49), (54), (55), (58), (59) yields 



(67) 
(68) 



(69) 



-F^ 77 (0) 



0.025 ±0.001 GeV -1 . 



F^ 77 (0) = 0.030 ± 0.001 GeV" 



(70) 



in excellent agreement with (69). Note that one can similarly compare the 7/, 7/ — > 77 
widths, as given by the relation (see, e.g., ref. [13]) 



r(p° - 77) 



F|o 77 (0)M 2 (P c 



64tt 



(71) 



with Fpo 77 (0) defined in (65), (67), (68), with those measured. Such a comparison 
gives (in keV): 0.51 ± 0.05 vs. |j 0.46 ± 0.04 for 77 -> 77, and 4.04 ± 0.27 vs. @] 
4.26 ± 0.19 for rj -> 77. 

Also, one can compare the values for the 77-77 mixing-independent i?-ratio, given 
by (using (65),(67),(68),(71)) 



R = 



r(rj -> 77) r(V -» 77) 



77 



7T" 



I J IT 1 n JI 



r(vr -> 77) 3 V /f 



/ 9 2 



With the values for f$/ f n and fg/ f n obtained above, our result is R = 2.25 ±0.15, vs. 
2.45 ± 0.35, as follows from using the experimentally measured masses and widths in 
the above expression. 
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5.2 ry, r/ — > tt + tt 7 decays 

These, as well as the P° — > 77, processes were extensively studied by Venugopal and 



Holstein [28fl in chiral perturbation theory. The analysis of experimental data for 



both of these processes done in ref. [|28| yields 

j- = 1.38 ±0.22, ^ = 1.06 ±0.03, 6 = (-22.0 ± 3.3)°, 

Jir Jit 

in good agreement with our Eqs. (59), (62), (63). 
5.3 J/ip — > 777, 77 7 decays 



These processes were studied by Kisselev and Petrov p7|. The values of fs/ fn 



f$l fir an d 9 extracted in ref. |E7| from the experimentally measured P° — > 77 and 



J j if) — > 777,777 widths, as given in the last column of Table I of ref. [27|, which 
corresponds to conventional mass-mixing angle relations, are 

j- = 1.12 ±0.14, A = 1.04 ±0.08, 9 = (-18.9 ±2.0)°, 

again in good agreement with our Eqs. (59), (62), (63). 

Thus, the three values of / 8 / f-x, hi fn and 9 agree with experiment (at least, as far 
as the processes considered above are concerned). As to the remaining f v / f n , fj/fm 
ratios also calculated in the paper, the experimental values of them, as extracted from 



data by the CELLO |34| and TPC/27 |35| collaborations, are, respectively, 

1.12 ±0.12, (72) 
1.06 ±0.10, (73) 



U 

fir 

k 



and 



k 



1.09 ±0.10, (74) 
0.93 ±0.09. (75) 



While the value calculated for f v / f n , Eq. (56), clearly agrees with both experimental 
values (72) and (74), the value calculated for f'/f^, Eq. (57), only marginally agrees 
with (73), and disagrees with (75) by almost 3 standard deviations. 

To clarify this point, let us note that the values of f v / f n and / / / f n were extracted 
by both CELLO and TPC/27 from experimental data on the transition form-factors 
T / /JO, — Q 2 ), assuming that the pole mass A / /\, which parametrizes their fits to 
the data, can be identified with 2ir\/2 f,> y Then, these pole fits to the data are 
presumed to join smoothly, as Q 2 — > 00, to the perturbative QCD predictions for 
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2^/ ^(0, — Q 2 ) PBfl , i.e., these fits would then agree with both the QCD asymptotic 
form ~ l/Q 2 and its coefficient. However, the values of both f v and / ' > quoted by 
the two groups, (in MeV) (94.0 ± 7.1, 89.1 ±4.9) [0, and (91.2±5.7, 77.8±4.9) |§, 
respectively, are all close to M(p)/(27ta/2) ~ 86.5 MeV, thus indicating on possible 
connection with the vector meson dominance interpretation of A^') ~ M(p) in the 
range of Q 2 investigated. 

On the other hand, as remarked in ref. | 57| , the model independent calculation 
by Gasser and Leutwyler p2| , testing the Goldberger-Treiman relations by Scadron 
p8| , and the calculation by Burden et al. P9fl , all agree that both and / ' should 
be noticeably larger than Our results f v ~ l-lf n , /„' ~ l-2/ w are in agreement 
with this. 

It therefore seems that the extraction of the values of f v and / / from the transition 
form-factors T^/^O, — Q 2 ) cannot be done accurately enough, at least in the range 
of Q 2 investigated so far. That this may indeed be the case is indicated by the 
experimental value f n a = 84.1 ± 2.8 MeV extracted by the same method, which 
is again close to M(p)/27r\/2. The central value of this f^o, 84.1 MeV, is ~ 10% below 
the well established value given in Eq. (49), 92.4 ±0.2 MeV. Such a large discrepancy 
cannot be explained by, e.g., small isospin violation, indicating therefore that the 
extracted values for both f v and / '< may well have been underestimated too. 



6 Concluding remarks 

As a lagniappe, we note that Eq. (50) may be combined with our extracted values 
for fg, f^, f i to obtain the value of A : 

A = 0.78 ±0.12 GeV 2 . 

This is consistent with the usual value 0.73 GeV 2 determined by the trace condition 



for (23), (25) without f s p 



Let us briefly summarize the findings of this work: 

i) We have found that many theoretical approaches suggest that the more natural 
object to study is the pseudoscalar mass squared matrix modified by the inclusion of 
the squared factors of the pseudoscalar decay constants. 

ii) We have shown that this modified mass squared matrix leads to new Schwinger's 
quartic mass and the Sakurai mass-mixing angle relations for the pseudoscalar meson 
nonet. 

iii) We have used these new relations for calculation of the pseudoscalar decay 
constants and mixing angle. We have demonstrated that, except where questions 
may be raised regarding the reliability of the extraction of the relevant quantities 
from direct experimental data, the results obtained are in excellent agreement with 
available data. 
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